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Presented  is  a  general  theoretical  framework  capable  of  describing  the  finite 
deformation  kinematics  of  several  classes  of  defects  prevalent  in  metallic  crystals. 
Our  treatment  relies  upon  powerful  tools  front  differential  geometry,  including 
linear  connections  and  covariant  differentiation,  torsion,  curvature  and 
anholonontic  spaces.  A  length  scale  dependent,  three-term  multiplicative 
decomposition  of  the  deformation  gradient  is  suggested,  with  terms  representing 
recoverable  elasticity,  residual  lattice  deformation  due  to  defect  fields,  and  plastic 
deformation  resulting  from  defect  fluxes.  Also  proposed  is  an  additional 
micromorphic  variable  representing  additional  degrees-of-freedom  associated 
with  rotational  lattice  defects  (i.e.  disclinations),  point  defects,  and  most 
generally,  Somigliana  dislocations.  We  illustrate  how  particular  implementations 
of  our  general  framework  encompass  notable  theories  front  the  literature  and 
classify  particular  versions  of  the  framework  via  geometric  terminology. 


1.  Motivation  and  background 

Methods  from  geometry  on  differentiable  manifolds  have  come  to  the  fore  in  recent 
years  for  describing  the  kinematics  of  elastoplasticity  [1-9].  An  understanding  of 
the  kinematics  of  defective  crystals  is  central  to  constructing  arguments  related  to 
nonlocal  or  gradient  theories  introducing  scale  effects  in  the  description.  As  such, 
the  primary  objective  of  the  present  work  is  to  provide  a  unified  framework  encap¬ 
sulating  the  kinematics  of  crystalline  incompatibilities  whose  descriptions  give  rise 
to  higher-order  gradient  continuum  plasticity  models.  Our  scope  is  restricted  to 
kinematics  of  defects;  for  complementary  treatments  of  material  forces,  the  reader 
is  referred  to  works  by  Eshelby  [10,  11],  Nabarro  [12],  Kroner  [13],  Maugin  [1-3] 
and  Gurtin  [14]. 

The  two-term  multiplicative  decomposition  of  the  deformation  gradient, 
F  =  FeFp,  was  proposed  by  Bilby  et  al.  [15]  and  Kroner  [16]  for  modeling  the 
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deformation  of  single  crystals  within  the  framework  of  continuously  distributed 
dislocations.  In  the  theory  of  Bilby  et  al.  [15],  F,  Fe,  and  Fp  are  described  as  the 
‘shape  deformation’,  the  ‘lattice  deformation’,  and  the  ‘dislocation  deformation’, 
respectively.  The  existence  of  a  generally  incompatible,  relaxed  intermediate 
configuration  was  described  in  the  years  prior,  for  example,  by  Eckart  [17], 
Kondo  [18],  Bilby  et  al.  [19]  and  Eshelby  [11].  Kondo  [18]  and  Bilby  et  al.  [19] 
are  usually  credited  with  discovering  the  relationship  between  the  density  of 
continuously  distributed  dislocations  and  Cartan’s  torsion  tensor  of  the  ‘crystal 
connection’,  which  we  will  see  may  be  defined  in  terms  of  a  spatial  gradient  of  the 
(inverse)  lattice  deformation.  The  geometry  induced  by  the  crystal  connection  is 
said  to  possess  a  ‘distant  parallelism’  or  a  ‘teleparallelism’  [20],  and  exhibits  the 
special  property  that  covariant  derivatives  taken  with  respect  to  this  connection  of 
spatially  constant  intermediate  configuration  slip  vectors  vanish.  The  null  curvature 
property  of  the  crystal  connection  (cf.  [21]),  together  with  knowledge  of  its  torsion 
and  either  the  covariant  plastic  or  elastic  strain  tensor  is  sufficient  to  uniquely  specify 
the  relaxed  intermediate  configuration,  as  demonstrated  by  Le  and  Stumpf  [4]. 
Le  and  Stumpf  [5]  and  Steinmann  [6]  also  provide  geometric  viewpoints  of  the 
kinematics  of  higher-gradient  finite  elastoplasticity  of  crystals  with  defects.  Also  of 
interest  is  the  description  of  Noll  [21]  admitting  multiple  reference  configurations. 
Rajagopal  and  Srinivasa  [22]  developed  a  framework  capable  of  characterizing 
evolving  elastic  symmetry  in  a  wide  class  of  fluids  and  solids,  wherein  the  deforma¬ 
tion  is  described  as  the  evolution  of  a  family  of ‘natural  configurations’.  Cermelli  and 
Gurtin  [23]  examined  the  kinematics  of  geometrically  necessary  dislocations  from  the 
perspective  of  multiplicative  elastoplasticity  with  plastic  deformation  characterized 
by  shearing  on  discrete  slip  systems.  Other  noteworthy  efforts  employing  differential 
geometry  to  describe  generalized  continua  include  the  manifold-based  treatment 
of  Marsden  et  al.  [24]  and  methods  founded  upon  gauge  theories  and/or  exterior 
calculus  [7,  25,  26]. 

Since  dislocations  associated  with  the  crystal  connection  are  required  to  sustain 
compatibility  of  the  total  deformation  and  may  induce  local  residual  stresses 
within  the  crystal,  they  are  often  labeled  ‘geometrically  necessary’  dislocations  in 
the  sense  of  Ashby  [27].  This  is  in  contrast  to  the  ‘absolute’  dislocations  [28,  29]  or 
‘statistically  stored’  dislocations  [27]  that  sustain  homogeneous  plastic  deformation. 
Lardner  [28]  showed  that  vanishing  of  the  torsion  tensor  of  the  crystal  connection 
corresponds  to  vanishing  of  excess  dislocations  of  the  same  sign.  Nye  [30]  associated 
a  dislocation  density  tensor  with  the  gradient  of  lattice  rotation,  i.e.  the  stress-free 
lattice  curvature.  Kroupa  [31]  proposed  a  loop  density  tensor  for  describing  the 
contributions  of  dislocation  loops  to  the  plastic  distortion  field  and  induced  internal 
stresses.  Kroner  [32,  33]  suggested  inclusion  of  various  statistical  moments  of 
dislocation  distributions  in  theories  of  work  hardening.  Grain  boundary,  surface, 
and  interface  dislocations  have  also  been  investigated  from  the  geometric  viewpoint 
[34-38], 

Kondo  [18,  39,  40]  proposed  more  general  frameworks  for  yielding  and  plastic 
deformation  allowing  nonzero  Riemann-Christoffel  curvature  associated  with  a 
connection  different  from  the  crystal  connection.  Kondo’s  frameworks  admit 
different  non-unique  anholonomic  (i.e.  incompatible)  intermediate  configurations, 
or  ‘tearings’  in  his  terminology,  for  the  same  crystal.  These  configurations  differ 
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from  the  unique  intermediate  configuration  favored  by  Bilby  et  al.  [19],  for  which 
the  curvature  of  the  crystal  connection  vanishes  and  teleparallelism  is  preserved. 
Kondo  [40]  associated  the  nonzero  curvature  with  the  presence  of  ‘rotational 
anomalies’  in  the  lattice,  which  could  be  interpreted  as  disclinations  [41,  42]. 
Connections  admitting  non-vanishing  curvature  tensors  have  since  been  used  by 
many  for  modeling  disclinations  in  Bravais  lattices  [43-48].  Rotational  defects  falling 
under  the  general  category  of  Somigliana  dislocations  (cf.  [10,  11,  49,  50])  were 
addressed  in  terms  of  a  path-dependent  integral  of  the  non-vanishing  curvature 
tensor  by  P^cherski  [51,  52].  Non-metric  connections  admitting  ‘extra -matter’  defects 
(e.g.  interstitial  atoms,  vacancies,  and  other  point  defects)  were  introduced  in 
continuum  theories  by  Kroner  [16,  53,  54],  Minagawa  [45]  and  De  Wit  [48]. 

The  remainder  of  the  present  paper  is  organized  as  follows.  Section  2 
defines  relevant  terminology,  followed  by  mathematical  preliminaries  in  section  3. 
In  section  4,  we  present  our  kinematic  description  at  the  macro-level,  specifically  the 
decomposition  of  the  deformation  gradient  for  a  local  volume  element  of  crystalline 
material.  Section  5  provides  our  kinematic  description  at  the  micro-level,  in  terms  of 
an  affine  connection,  possibly  non-metric  with  a  potentially  non-vanishing  curvature 
tensor.  We  illustrate  throughout  sections  4  and  5  how  our  description  encompasses 
the  physics  of  translational  dislocations,  disclinations,  point  defects,  and  Somigliana 
dislocations,  and  we  discuss  how  particular  versions  of  our  framework  reduce 
to  noteworthy  models  in  prior  literature.  Key  features  are  summarized  in  tabular 
form  in  section  6. 

2.  Terminology 

Here,  we  summarize  terminology  that  is  often  a  source  of  confusion  to  readers 
not  well-versed  in  differential  geometry  and  its  application  to  continuum  mechanics. 
The  definitions  given  immediately  below  are  free  of  notation  and  hence  are 
somewhat  qualitative;  in  many  cases  corresponding  equations  may  be  found 
later  in  the  paper,  or  the  reader  may  elect  to  consult  more  extensive  sources 
(e.g.  [55-57]). 

•  Configuration.  A  configuration  is  defined  as  a  generally  time-dependent 
realization  of  a  body.  In  finite  deformation  theories,  one  often  speaks  of  the 
‘reference’,  ‘initial’,  ‘undeformed’,  or  ‘Lagrangian’  configuration,  usually  and 
somewhat  arbitrarily  taken  at  zero  time.  Similarly,  the  ‘spatial’,  ‘current’, 
‘deformed’,  or  ‘Eulerian’  configuration  is  usually  taken  at  the  current  time. 
These  two  configurations  are  assumed  to  be  holonomic  to  one  another, 
meaning  that  current  coordinates  can  be  written  as  smooth,  single-valued 
functions  of  reference  coordinates,  and  vice-versa.  In  contrast,  the  virtual 
‘intermediate’,  ‘relaxed’,  or  ‘unloaded’  configuration  of  elastoplasticity  is 
often  anholonomic,  since  its  coordinates  cannot  always  be  prescribed  as 
single-valued  functions  of  reference  or  current  coordinates.  The  reader  is 
referred  to  Wang  and  Truesdell  [56]  for  a  more  in-depth  discussion. 

•  Metric  tensor.  A  metric  tensor  (or  simply  a  ‘metric’)  is  a  rank  2  covariant 
tensor  that  defines  the  scalar  product  of  contravariant  vectors,  and  conse¬ 
quently,  the  squared  length  of  a  vector.  A  metric  tensor  that  is  both  symmetric 
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and  positive  definite  is  called  a  Riemannian  metric  tensor f;  non-Riemannian 
metrics  will  not  be  employed  here.  Metric  tensors  perform  other  consequential 
functions,  such  as  raising  and  lowering  indices  and  defining  scalar  products. 

•  Connection.  A  linear:}  connection  is  a  rank  three  construct  that  defines 
a  covariant  derivative  of  vectors  (cf.  (1))  and  tensors  of  higher  rank  (4). 
Alternatively  and  equivalently,  the  covariant  derivative  operation  defines  the 
connection  coefficients  (i.e.  Christoffel  symbols).  See  Eisenhart  [58]  or  Boothby 
[59]  for  additional  perspective. 

Note  that  a  particular  configuration  or  space  can  be  assigned  more  than 
one  connection,  just  as  it  can  be  assigned  more  than  one  metric  tensor.  We  can 
often  classify  the  pair  of  {configuration,  connection}  or  triad  of  {configuration, 
connection,  metric  tensor}  as  one  or  more  of  the  following  five  types  of 
‘geometric  spaces’: 

(1 )  Euclidean  space.  For  a  space  to  be  classified  as  Euclidean  or  non- Euclidean,  it 
must  include  a  configuration,  a  metric  tensor,  and  a  connection.  A  Euclidean 
space  satisfies  three  requirements:  (i)  the  torsion  tensor  from  the  connection 
vanishes,  (ii)  the  covariant  derivative  of  the  metric  tensor  vanishes,  and 
(iii)  the  Riemann-Christoffel  curvature  tensor  from  the  connection 
coefficients  vanishes.  A  Euclidean  //-space  permits  at  each  location  a 
transformation  from  local  coordinates  to  a  global  /(-dimensional  Cartesian 
coordinate  system.  In  continuum  mechanics,  both  the  reference  and  current 
configurations  are  typically  viewed  as  three-dimensional  Euclidean  spaces, 
with  the  motion  acting  as  a  diffeomorphism  (i.e.  a  differentiable  homeomor- 
phism ,  or  a  differentiable  one-to-one,  invertible  mapping)  between  these  two 
configurations  [61].  Notice  the  distinction  between  Euclidean  and  Cartesian-. 
the  former  (e.g.  spherical  coordinates  in  three  dimensions)  admits  a  global 
transformation  to  the  latter  (e.g.  three  constant  and  orthonormal  basis 
vectors). 

(2)  Anholonomic  space.  An  anholonomic  space  is  regarded  as  a  configuration 
defined  by  a  non-integrable,  two-point  deformation  map.  In  multiplicative 
plasticity,  since  Fe_1  and  Fp  are  generally  non-integrable  or  ‘anholonomic’ 
functions  of  current  and  reference  coordinates,  respectively,  the  correspond¬ 
ing  intermediate  configuration  is  generally  anholonomic.  Coordinate 
functions  on  such  an  anholonomic  space  are  necessarily  discontinuous, 
multi-valued  functions  of  current  or  reference  configuration  ‘holonomic’ 
coordinates.  Anholonomicity  is  related  to  Cartan’s  torsion  tensor  of  a 
certain  connection  (i.e.  the  crystal  connection)  constructed  from  the  lattice 
deformation  in  our  theory,  as  discussed  in  section  5.1  (see  equation  (21)). 
Anholonomic  coordinates  are  discussed  at  length  by  Schouten  [55]. 


f  A  non-Riemannian  metric  tensor  would  not  be  positive-definite.  Such  metrics  arise  in  more 
generalized  geometries  such  as  Finlser  spaces,  an  example  of  which  is  Minkowski’s  spacetime 
(cf  Rund  [60]).  Metrics  are  symmetric  by  definition. 

}  In  the  present  work  we  focus  on  linear  (also  called  affine)  connections,  which  strictly 
obey  (1).  Nonlinear  connections  arise  in  more  generalized  spaces  such  as  Finsler 
manifolds  [60]  and  will  not  be  dealt  with  here. 
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(3)  Cartan  space.  A  configuration  with  a  connection  admitting  a  non-vanishing 
Cartan’s  torsion  tensor  is  labeled  a  ‘Cartan’  space,  sometimes  called  a 
‘non-symmetric’  space.  We  need  only  consider  the  pair  {configuration, 
connection}  to  label  a  space  as  Cartan  or  non-Cartan. 

(4)  Riemannian  space.  For  a  space  to  be  labeled  as  ‘Riemannian’  or  ‘non- 
Riemannian’  it  must  include  the  pair  {configuration,  connection}  —  a  metric 
is  not  needed  for  our  designation.  A  Riemannian  space  is  defined  here  as  a 
configuration  with  connection  coefficients  whose  components  yield  a  nonzero 
Riemann-Christoffel  curvature  tensor.  A  space  with  a  non-vanishing  curva¬ 
ture  is  non-Euclidean:  for  example,  a  global  2D  Cartesian  coordinate  system 
cannot  be  used  to  parameterize  a  shell  unless  the  shell  is  flat.  Conversely, 
a  space  with  vanishing  curvature  is  ‘non-Riemannian’,  and  the  connection 
is  said  to  be  ‘integrable’. 

(5)  Non-metric  space.  For  one  to  label  a  generalized  space  as  ‘metric’  or  ‘non¬ 
metric’  the  configuration  must  be  assigned  both  a  linear  connection  and  a 
metric  tensor;  we  must  examine  the  trio  of  {configuration,  connection,  metric 
tensor}.  In  a  non-metric  space,  the  covariant  derivative  of  the  metric  tensor 
taken  with  respect  to  the  connection  is  nonzero. 

The  above  terminology  is  not  fully  consistent  in  the  literature  [6,  19,  21,  40,  55,  57], 
We  have  attempted  to  use  terminology  appearing  most  frequently  and  logically  for 
describing  crystals. 


3.  Mathematical  preliminaries 

Reviewed  next  are  fundamental  geometric  concepts.  Definitions  are  given  first, 
followed  by  a  discussion  of  strain  compatibility  and  anholonomic  coordinates. 
For  more  in-depth  treatments,  the  reader  is  referred  Schouten  [55],  Wang  and 
Truesdell  [56]  and  Marsden  and  Flughes  [57],  Boldface  type  is  used  for  vectors 
and  tensors  of  higher  rank,  while  scalars  and  components  of  tensors  are  written 
in  italics.  In  order  to  appeal  to  a  broad  audience  (e.g.  engineers,  physicists,  and 
mathematicians)  we  often  use  the  index  notation  favored  by  earlier  differential 
geometers  [55,  58]  even  though  many  expressions  could  be  presented  more  concisely 
via  more  modern  direct  notation  of  tensor  analysis  on  manifolds  and  exterior 
calculus  [7,  57]. 

3.1.  Definitions 

An  n-manifold  is  a  set  M  such  that  for  each  P  e  M  there  is  a  subset  U  of  M  contain¬ 
ing  P,  and  a  one-to-one  mapping  from  U  onto  an  open  set  in  D?"  (the  n-  space  of 
real  numbers).  Changes  in  coordinates  over  coincident  regions  are  assumed  to  be 
infinitely  differentiable  if  the  manifold  is  smooth.  For  the  present  discussion,  we 
introduce  a  smooth  manifold  B0 ,  which  we  associate  with  the  reference  configuration 
of  a  body  in  the  terminology  of  finite  deformation  continuum  mechanics.  Manifold 
B()  is  parameterized  by  coordinates  XA  (A  —  1,2,3),  with  the  tangent  space  at 
each  X  spanned  by  the  natural  basis  vectors  =  dAX,  where  dA  =  d/dXA. 
The  ensemble  of  local  tangent  spaces  of  B(),  or  tangent  bundle ,  is  written  TB0. 
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The  metric  tensor  corresponding  to  natural  basis  vectors  on  B0  is  denoted  G,  where 
Gab  —  Ga-Gb. 

A  linear  connection  (also  often  called  an  affine  connection)  on  a  manifold  BQ 
induces  an  operation  V  that  assigns  to  two  vector  fields  V,  W  e  TB0  a  third  vector 
field  VvW  e  TB0,  called  the  covariant  derivative  of  W  along  V,  such  that 

(i)  VvW  is  linear  in  both  W  and  V; 

(ii)  V/vW  =  /VvW  for  scalar  functions  /;  (1) 

(hi)  Vy(/W)  =  /V  y  W  +  V  (D/)W . 

In  (iii)  above,  V(D/)  =  VAd Af  The  Christoffel  symbols  (or  connection  coefficients) 
rBc  °f  a  connection  V  on  B0  are  defined  on  a  coordinate  system  XA,  with  holonomic 
basis  { G_4 } ,  by 

FbcGa  —  VGbGc.  (2) 

Using  equations  (1)  and  (2),  the  covariant  derivative  of  W  along  V  is  written  as 

Vv  W  =  ( VBdB  WA  +  Pbac  Vb  U^)  Ga .  (3) 

An  affine  connection  on  the  tangent  bundle  TB()  defines  parallel  transport  of  vectors 
on  different  tangent  spaces  in  TB0,  such  that  vectors  in  different  tangent  spaces  may 
be  compared,  permitting  calculation  of  the  change  of  constituent  vectors  of  a  vector 
field  over  the  entire  tangent  bundle.  A  vector  is  said  to  undergo  parallel  transport 
(with  respect  to  a  connection  with  covariant  derivative  V)  along  paths  for  which  its 
covariant  derivative  vanishes.  For  example,  a  vector  W  is  considered  to  be  “parallel” 
along  a  curve  k(t)  if  VvW  =  0,  where  V  =  31/3 1  is  the  tangent  vector  to  the  curve 
at  the  current  value  of  parameter  t. 

The  covariant  derivative  is  applied  to  tensor  fields  as  [55]: 

V7  A(V\A-F  o  aA...F  I  y^.-A  aRB...F  I  I  j~'--F  aA...ER 

v Na\a)g...M  —  (jN/iG...M  ~r  1  NRAG...M  “t  +  1  NR^G-.M 

(4) 

r..R  jA...F  r..R  .  •(.../•' 

—  1  NGaRH...M  ~  —  1  NMaG...LR > 

where  we  have  placed  the  index  of  covariant  differentiation  as  a  subscript 
immediately  following  the  V-operator.  Analogous  to  parallel  transport  of  vectors, 
a  tensor  is  parallel  transported  along  a  curve  if  its  covariant  derivative  vanishes  as 
the  tensor  is  ‘dragged’  along  the  curve. 

The  torsion  tensor  T  of  a  connection  is  defined  by 


T(V,  W)  =  VvW  —  VWV  —  [V,  W],  (5) 

where  the  Lie  bracket  of  vector  fields  V  and  W  on  TB0  is  given  by 

[V,  W]  =  ( VBdBWA  -  WBdBVA) Ga  (6) 

Using  Equations  (3),  (5)  and  (6),  the  coordinate  representation  of  the  torsion  tensor 
is  found  as 


T  =  {r-Bc  -  rd) Ga  ®  Gfi  ®  Gc  =  2r[fAc]GA  ®  G5  ®  Gc, 


(7) 
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where  indices  enclosed  in  brackets  [•]  are  anti-symmetrized,  i.e.  2AyAB\  —  Aab  —  Aba, 
and  ®  is  the  outer  product.  The  torsion  tensor  of  a  linear  connection  is  often  called 
Cartan’s  torsion ,  by  association  with  geometer  E.  Cartan  [62].  A  connection  with 
vanishing  torsion  is  said  to  be  symmetric.  For  a  smooth  manifold  Bq  with 
Riemannian  metric  G,  there  is  a  unique  affine  connection,  which  we  label  now 
on  Bq  that  is  torsion-free  (p  —  o)  and  for  which  parallel  transport  preserves  dot 
products  of  vectors,  i.e.  y  G  =  0-  It  is  called  the  Levi-Civita  connection  [57]  or  the 
Riemannian  connection  [63],  The  connection  coefficients  of  the  covariant  derivative  y 
are  defined  by 


T  =  -  GAB(dBGcB  +  9cGdb  ~  9bGBc)Ga  <®  Gc  ®  G°.  (8) 

A  connection  yielding  a  null  covariant  derivative  of  the  metric  tensor  is  called  a 
metric  connection.  The  Levi-Civita  connection  of  (8)  is  metric.  Note  that  in  Cartesian 
coordinates,  GAB  =  8AB  =  constant,  meaning  that  the  Christoffiel  symbols  from  GAB 
vanish  by  (8). 

The  Riemann—Christojfel  curvature  tensor  associated  with  an  arbitrary  linear 
connection  with  covariant  derivative  V  on  B0  is  a  {3}  tensor  satisfying  [57] 

R(a,  W,,  W2,W3)  =  a(VW|VW2W3  -  VW2VW|W3  -  V[w„w2]W3),  (9) 

where  a  and  W,-  denote,  respectively,  arbitrary  one-forms  (i.e.  covariant  vectors)  and 
contravariant  vectors  in  configuration  B0.  The  component  representation  of  R  is 
given  by 


R  =  ( dcr-64B  -  dDrci  +  rcAErDEB  -  rDAErfB)GA  ®  gb  ®  gc  ®  g".  (io) 

G 

We  denote  by  R  the  curvature  tensor  formed  by  inserting  the  Levi-Civita  con¬ 
nection  (8)  into  (10).  A  geometric  space  BQ  with  metric  G  having  r  =  0  is  called  flat. 
One  may  show  (cf  [55])  that  r  —  0  if  and  only  if  we  may  assign  coordinate 
bases  at  each  X  e  B0  such  that  GAB  $ab',  i-e-  if  and  only  if  B0  is  Euclidean.  In  fact, 
r  —  0  are  compatibility  conditions  for  existence  of  connection  coefficients  p  derived 
from  a  Euclidean  metric  G  via  (8). 


3.2.  Finite  strain  compatibility 

Assume  that  the  coordinates  XA  covering  B0  undergo  a  smooth,  time-dependent 
mapping  to  ‘current’  coordinates  xa  —  tpa(XA,i)  spanning  B.  Current  configuration 
B  is  parameterized  by  spatial  coordinates  xa  (a  =  1,2, 3),  with  the  tangent  space  at 
each  point  on  B  spanned  by  the  natural  basis  vectors  ga  =  3flx.  The  tangent  bundle  in 
the  current  configuration  is  TB.  and  the  metric  tensor  g  corresponding  to  the  natural 
bases  on  B  has  components  gah  —  gfl-g b.  The  linear  tangent  map  F  =  Tip  is  labeled 
the  local  deformation  gradient'. 

F=T“g„®G"  =  f^ga®G". 


(11) 
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The  right  Cauchy-Green  deformation  tensor  is  the  pull-back  of  the  current 
configuration  metric: 

C  -  <p*(g)  =  FaAgahFhBGA  ®  Gb,  (12) 

where  <p*  denotes  the  pull-back  operation  by  F  (cf.  [57]).  The  strain  tensor  C  assigns 
any  (infinitesimal)  vector  dX  —  dXAGA  the  length  it  obtains  after  deformation 
by  F,  i.e.  ||t/x||  =  ||FrfX||  =  V(C dX,  dX)  =  V" CABdXAdXB.  We  introduce  a  linear 
connection 


r  =  -(C  ')  (dDCCB  +  dcCDB  —  dBCDC)GA  ®  Gc  ®  G° .  (13) 

c  c 

We  then  form  the  curvature  R  by  substituting  the  components  of  T  into  (10).  Since 
C  =  <p*(g),  one  may  show  [57],  using  the  propertiesc,of  a  linear  connection  (3.1)  and 
the  definition  of  the  curvature  tensor  (3.9),  that  R  =  <p*(R).  Hence,  if  the  tensor 
field  C  is  derivable  from  a  motion  function  <p(X,  t)  and  if  R  =  0,  then  C  is  compatible 
and  r  —  tp* (0)  =  0-  Notice  that  C-compatibility  does  not  require  specification  of  a 
unique  current  configuration,  since  C  is  independent  of  the  rotation  tensor  Q 
associated  with  the  right  polar  decomposition  F  =  QU.  Furthermore,  as  discussed 
by  Truesdell  and  Noll  [101]  and  Acharya  and  Bassani  [100],  even  if  C  is  compatible, 
the  integrability  of  an  arbitrary  deformation  gradient  field  F  generating  the 
strain  field  C  may  be  destroyed  by  rotations  Q  that  do  not  arise  from  rigid-body 
transformations  of  <p.  A  material  in  such  a  condition  is  said  to  be  in  a  state  of 
contorted  aleotropy  in  the  sense  of  Noll  [21].  The  converse  of  the  previous  theorem 
has  also  been  proven,  albeit  only  locally  [58].  In  other  words,  given  a  positive- 
definite,  symmetric,  second-rank  tensor  C  whose  curvature  vanishes  (i.e.  r  —  o) 
then  at  any  point  with  coordinates  XA  e  B0  there  exists  a  neighborhood  U0  of  XA 
endowed  with  a  mapping  <p :  B0  D  U0  ->  U  C  B  whose  deformation  tensor  is  C.  One 
can  consider  C-compatibility  an  outcome  of  deformation  gradient  compatibility 
(F-compatibility),  which  is  discussed  next. 


3.3.  Deformation  gradient  compatibility 

Consider  now  a  field  of  contravariant  (basis)  vectors  {§„},  where  (a—  1,2,3), 
spanning  configuration  space  B.  By  introducing  the  two-point  deformation  map 
F  =  F°  A  ga  ®  G4,  which  we  assume  is  defined  globally,  we  can  push  forward  vectors 
V  e  TBq  to  the  tangent  bundle  TB: 

FV  =  F*Aga  ®  G4  (  VbGb )  =  PU  VBgu{GA,  GB)  =  FA  VBSjga  =  F*A  VAga  e  TB.  (14) 


The  basis  vectors  are  tangent  to  globally-continuous  coordinate  curves  xa  (i.e. 
gff  =  dax  for  some  coordinate  parameterization  x“(X))  if  and  only  if  the  following 
integrability  conditions  hold  for  F  [55]: 


—  9aF"b 


d2xa 

dXAdXB 


d2xa 

dXBdXA  ' 


(15) 
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If  the  conditions  (15)  are  not  satisfied,  then  {§„}  is  called  an  anholonomic  coordinate 
basis,  the  deformation  map  F  is  called  an  incompatible  map,  and  the  configuration 
space  B  is  called  an  incompatible  configuration  or  an  anholonomic  space.  Since 
conventional  partial  differentiation  with  respect  to  anholonomic  coordinates  does 
not  apply,  we  define  differentiation  with  respect  to  the  xa  as 


9(0 

9x“ 


p~lA 
dxA  •“ 


9«(‘)  =  dA(-)F~1A. 


(16) 


A  Levi-Civita  connection  (i.e.  a  connection  both  torsion-free  and  metric  with 
respect  to  g,  the  latter  with  components  gap  =  ga  •  gf  on  B  may  not  exist,  since 
the  field  of  tangent  vectors  ]gff]  may  not  be  sufficiently  smooth  over  all  of  B  to 
admit  coordinate  differentiation  (with  respect  to  coordinates  that  may  in  fact  not 
exist,  i.e.  anholonomic  coordinates).  However,  in  the  anholonomic  (intermediate  or 
natural)  configurations  of  elastoplasticity  theory,  each  local  volume  element  is  often 
referred  to  an  external  system  of  coordinates  with  Euclidean  metric  tensor,  typically 
taken  as  a  global  Cartesian  frame  [64,  65],  although  such  an  assumption  is  not 
necessary  from  a  physical  perspective  [2,  66].  The  Riemann-Christoffel  curvature 
formed  from  the  covariant  strain  measure  CAB  =  FaAgapFpB  does  not  necessarily 
vanish  (unless  holonomic  xa  are  available),  in  contrast  to  the  curvature  tensor 
derived  from  connection  (13)  formed  from  the  deformation  C. 


4.  Deformation  gradient  kinematics 

Our  framework  encompasses  kinematic  descriptions  at  two  length  scales,  which  we 
shall  refer  to  henceforth  as  the  macroscale  and  the  microscale.  The  macroscopic 
description,  predicated  upon  the  multiplicative  decomposition  of  the  deformation 
gradient,  describes  the  elastic-plastic  kinematics  for  a  volume  element  of  crystalline 
material  in  an  average  sense.  The  microscopic  description,  predicated  upon  an 
additively-decomposed  linear  connection,  embodies  higher  spatial  gradients  of 
lattice  deformation  as  discussed  in  section  5,  enabling  one  to  interpolate  for  lattice 
stretch  and  orientation  between  centroids  of  neighboring  elements. 

Let  F  denote  the  volume-averaged  deformation  gradient  (equation  (11))  for  an 
element  of  crystalline  material.  Most  generally,  we  decompose  F  multiplicatively  as 

F  =  Tip  =  —  =  FiFp  =  FeF'Fp,  (17) 

9X 

where  V1  =  FeF‘  represents  the  total  lattice  deformation,  Fe  is  the  recoverable  elastic 
deformation,  and  Fp  is  the  plastic  deformation,  leaving  F1  as  the  residual  lattice 
deformation  associated  with  defects  and  other  sources  of  heterogeneity.  None  of 
the  mappings  Fe,  F1,  Fp,  or  VL  need  be  holonomic  when  considered  from  the  scale 
of  neighboring  volume  elements.  The  inverse  of  the  recoverable  elastic  deformation, 
Fe_1,  is  achieved  via  instantaneous  traction  removal  from  the  surface  of  the  crystal¬ 
line  element,  and  may  include  rigid  body  rotation  of  the  solid  as  well  as  recoverable 
thermal  deformation  (e.g.,  average  contraction  of  the  lattice  upon  cooling  to 
a  reference  temperature).  Since  lattice  defects  are  not  introduced  during  this 
hypothetical  unloading  process,  Fe  does  not  alter  the  holonomicity  (or  lack  thereof) 
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of  the  material  within  the  volume  element.  However,  since  the  deformation 
achieved  upon  unloading  may  not  be  compatible  across  neighboring  elements,  Fe 
is  anholonomic  at  the  macroscale.  The  deformation  Fp  embodies  contributions  of 
defects  that  leave  the  lattice  unperturbed  in  an  average  sense,  e.g.  dislocation  glide. 
Finally,  F1  accounts  for  remaining  physics,  including 

•  Residual  strain  and  rotation  fields  attributed  to  non-redundant  (i.e.  geometri¬ 
cally  necessary)  dislocations  [67] 

•  Rotation  of  the  lattice  attributed  to  rotational  defects  (i.e.  disclinations), 
which  may  describe  subgranular  defect  walls  formed  during  large  deforma¬ 
tions  (cf.  [68])  or  deformation  twins  [69,  70] 

•  Residual  thermal  strains  due  to  heterogeneous  thermal  expansion  properties 
within  the  volume  element,  perhaps  induced  by  defects  in  the  microstructure 

•  Volumetric  defects  such  as  interstitials  or  voids  resulting  in  lattice  expansion 
or  contraction 

•  Phase  transformations  that  alter  the  lattice  arrangement  [70] 

Decomposition  (17)  and  the  F'  deformation  mapping  were  introduced  by  the 
present  authors  to  describe  elastoplastic  polycrystals  [71]  and  elastoplastic  single 
crystals  [67],  both  in  the  context  of  dislocation-based  crystal  plasticity  theory. 
The  concept  is  generalized  here  to  account  for  various  residual  lattice  deformation 
modes  listed  above.  Note  that  (17)  is  scale  dependent,  depending  upon  the  size  of 
the  crystal  volume  element  to  which  it  is  applied,  and  resolution  dependent,  depend¬ 
ing  upon  the  way  in  which  local  kinematics  within  the  volume  element  are  resolved 
and  accounted  for  in  the  description.  For  example,  the  relative  magnitudes  of 
components  of  F1  will  differ  for  a  single  crystal  containing  a  few  non-redundant 
dislocations  and  a  heterogeneously-deforming  polycrystal  supporting  large  residual 
stresses  due  to  intergranular  incompatibility.  Kratochvil  [72]  postulated  a  three-term 
decomposition  similar  to  (17),  although  the  length  scale  dependence  was  not 
emphasized.  Our  theory  reduces  to  one  suggested  by  Bilby  and  Smith  [99]  for 
residual  elastoplasticity  when  Fe  =  1  and  the  elastic  unloading  corresponds  to  global 
unloading,  leading  to  the  reduced  decomposition  F  =  F'FP  (see  [73]).  And  our 
decomposition  reduces  to  that  of  Lardner  [44,  74]  and  Pe$herski  [51,  52]  when  we 
set  F'  =  R,  where  R  is  an  orthogonal  matrix  describing  lattice  rotation  caused  by 
disclinations. 

Macroscopic  configurations  of  a  crystalline  element  corresponding  to  (17)  are 
depicted  in  figure  1 .  Intermediate  configurations  B  and  B  are  generally  anholonomic, 
spanned  by  anholonomic  coordinates  xa  and  xa,  respectively,  while  reference  and 
spatial  configurations  B0  and  B  are  holonomic,  spanned  by  continuous  coordinates 
XA  and  x“ .  respectively.  Configuration  B  supports  nonzero  traction  vector  t,  while 
configurations  B  and  B  are  traction  free,  i.e.  t  =  0  and  t  =  0.  Lattice  director  vectors 
are  mapped  across  configurations  via 

da  =  (FL~ ')%,  d„  =  (F-')%  da  =  (18) 

while  Fp  is  assumed  to  leave  the  lattice  unperturbed.  When  the  lattice  directors 
are  heterogeneous  within  the  crystalline  volume  element  (e.g.  in  polycrystals  or  in 
subdivided  single  crystals),  the  triads  of  (18)  and  figure  1  are  understood  to 
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Figure  1.  Deformation  mappings  and  configurations  for  crystalline  volume  element. 

represent  suitably-defined  ‘average’  directors  for  the  volume  element  [75].  In 
non-cubic  lattices,  the  directors  of  (18)  are  not  necessarily  parallel  to  primitive 
translation  vectors  of  the  unit  cell,  but  are  instead  an  external  construct. 
Upon  assuming  that  the  directors  in  B  are  orthonormal,  and  assigning 
an  external  Cartesian  coordinate  system  to  anholonomic  space  B  with  metric 
tensor 


Aa-Aft  =  %afS  =  Saf  (19) 

we  define  a  metric  tensor  CL  referred  to  configuration  B  and  associated  with  the 
director  strain: 


Cat,  =  d„  •  d/,  =  FLn~laFL-lf>  d, 


j_  c’L— loto 

a  u p  —  r.a  °af)r .b 


(20) 


The  metric  CL  will  play  a  prominent  role  in  the  continuation  of  our  framework 
in  section  5. 

The  particular  sequence  of  deformation  maps  in  (17)  is  chosen  based  on 
rational  arguments.  The  plastic  deformation  term  Fp  is  the  rightmost  in  our 
decomposition,  as  it  leaves  the  lattice  vectors  unaltered,  permitting  identification 
of  configuration  b  of  figure  1  with  Mandel’s  [76]  isoclinic  configuration.  The 
residual  deformation  due  to  micro-heterogeneity  in  the  presence  of  lattice  defects, 
F1,  is  placed  second,  as  it  affects  the  lattice  directors  via  (18),  yet  is  assumed 
unaffected  by  superposed  rigid  body  motion  or  a  change  in  spatial  coordinate 
frame  [72].  Finally,  the  recoverable  elastic  deformation  Fe  =  VeRe  logically 
assumes  the  leftmost  position,  as  the  stretch  Ve  is  associated  with  unloading 
of  the  average  traction  acting  on  the  volume  element  from  the  current  state 
and  the  orthogonal  matrix  Re  accounts  for  all  lattice  rotations  not  embodied 
in  F1,  including  net  rigid  body  motions  of  the  solid.  The  rotational  components 
R'  of  F'  and  Rp  of  Fp  are  assumed  to  evolve  independently  of  rigid  body  motions 
of  the  solid.  Although  numerous  decompositions  are  possible  ( cj .  [77]),  we  find 
our  description  (17)  most  appropriate  and  realistic  from  the  physical  standpoint. 
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Figure  2.  Configurations  of  crystal  element  containing  edge  dislocation:  (a)  Volterra  model 
and  (b)  lattice  model. 


Figure  2  describes  the  physics  of  equation  (17)  from  the  standpoint  of  (a)  a 
Volterra  process  and  (b)  a  cubic  crystalline  lattice  for  a  volume  element  contain¬ 
ing  a  single  edge  dislocation.  The  tangent  line  of  the  defect  is  written  as  £0,  and 
the  Burgers  vector  is  labeled  b,  with  adornments  corresponding  to  the  particular 
configuration.  Notice  how  the  displacement  discontinuity  in  the  wake  of  the 
dislocation  is  sealed  by  the  residual  lattice  deformation  F'.  Prior  to  the  applica¬ 
tion  of  F1,  for  the  body  labeled  B  in  figure  2,  elastic  strain  fields  are  absent, 
and  the  plastic  deformation  gradient  Fp  is  completely  defined  in  terms  of  the 
relative  motion  of  the  two  halves  of  the  lattice  on  opposite  sides  of  the  slip 
plane  (dotted  line  in  figure  (2b)).  Note  also  that  B  is  generally  locally  incom¬ 
patible  across  the  entire  slip  path,  contains  slip  discontinuities  but  no  lattice 
strains,  and  is  free  of  internal  residual  stress  fields.  This  is  in  contrast  to  5, 
which  may  contain  internal  residual  stress  and  lattice  strain  fields  but  includes 
no  discontinuities  except  for  those  in  the  immediate  vicinity  of  individual  defect 
lines  contained  within  the  volume. 


5.  Micromorphic  kinematics 

At  the  microscopic  level,  spatial  changes  in  the  lattice  directors  d„  are  described  by 
parallel  transport  with  respect  to  a  linear  connection.  The  absolute  change  of  the 
director  vector  field  in  the  current  configuration  is  found  in  terms  of  the  covariant 
derivative  V,  i.e. 


—  da  b  —  nX,  (21) 

with  the  subscripted  comma  denoting  partial  coordinate  differentiation  with  respect 
to  x.  The  Christoffel  symbols  of  the  connection  are  defined  by  [45,  46] 

_  irLa  j-'L—la  ,  /~ \..a  7 -<..a  < 

1  cb  =  r  a  r.b,c  +  \lcb  ~  1  cb  \lcb^ 


(22) 
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where  rcb  =  F^'F^f101  —  —F^cF^~la  are  coefficients  of  the  crystal  connection  of 
non-Riemannian  dislocation  theories  [4,  5,  19,  21],  and  Q'cb  are  micromorphic 
degrees-of-freedom  representing  contributions  of  distributed  defects  to  the  spatial 
gradients  of  the  lattice  director  field.  Upon  assuming  that  Vftda  =  0  [46],  the  connec¬ 
tion  (22)  allows  one  to  interpolate  for  the  directions  and  magnitudes  of  the  lattice 
directors  between  centroids  of  neighboring  crystal  volume  elements.  For  example,  in 
the  trivial  situation  when  F'cb  =  0,  the  lattice  directors  are  spatially  constant.  The 
crystal  connection  (T)  component  of  (22)  accounts  for  effects  of  first-order  spatial 
gradients  of  (the  inverse  of)  F  ,  while  the  micromorphic  variable  Q  of  (22)  accounts 
for  additional  spatial  variations  of  lattice  directors  not  captured  by  the  first  gradient 
of  F  .  The  average  continuum  deformation  of  the  director  vectors  located  at  the 
volume  element’s  centroid  is  determined  by  F  ,  as  indicated  by  equation  (18).  Lattice 
stretch  and  rotation  gradients  (e.g.  at  the  scale  of  subgrain  cells  and  cell  blocks 
( cf.  [68,  78])  are  represented  by  the  coefficients  of  equation  (22). 

The  covariant  components  of  Q  are  defined  via  lowering  by  the  metric  CL 
of  equation  (20): 

Qcba  -  QcbCLda.  (23) 

Quantities  Qcba  and  Q'b  are  effectively  equivalent  micro-deformation  measures 
only  when  lattice  strains  are  negligible,  i.e.  when  Cda  8da.  The  covariant  derivative 
of  CL  is  found  as 

VcCLab  =  C£c  -  f£aCLdh  -  t£cLad -Q£cLdb  -  Q£cLad  =  -2 Qc(ah),  (24) 

' 

=0 


where  parentheses  denote  symmetrization,  e.g.  2Qc^hd)  =  Qcba  +  Qcab.  Components 
of  the  torsion  tensor  T  of  the  connection  T  are  given  by  (see  equation  (7)) 

hi  -  hi  -  hi  =  hi + 2Qm  *  (25) 

where  Tch  is  hie  torsion  of  the  crystal  connection.  The  components  of  the 
Riemann-Christoffel  curvature  tensor  Rh‘!d  formed  from  the  connection  coefficients 
rcb  are  found  using  equation  (10): 

n-  -a  _ 

Kbcd  =  1  db,  c  ~  1  cb,  d  '  1  ce  1  db  ~  1  de 1  cb> 

which,  because  t'cb  =  t'cb  +  Q'cb,  we  are  able  to  rewrite  as  [55] 

hcad  -  h£  +2 %Q-&  +  QclQfb  -  Q'deQ'cb  +  hdQel  (27) 

=0 


where  the  curvature  from  the  crystal  connection,  Rbfd,  vanishes  identically  as  shown, 
since  f  is  integrable  (see  equation  (53)).  The  fully  covariant  curvature 

hhcd  =  Cafhcd  *S  t^5] 

R[ab]cd  —  IFJycQd^bal  +  F'cdQe[ba ],  R{ab)cd  —  2V \cQ^(ab)  +  TfdQe{aby  (28) 
We  see  from  (28)  that  R  vanishes  completely  when  Q  =  0. 


3996 


J.  D.  Clayton  et  al. 


The  concept  of  a  Somigliana  dislocation  [79]  can  be  used  to  characterize  the 
displacement  discontinuities  from  translational  and  rotational  dislocations  (i.e. 
disclinations)  of  Volterra  [41],  point  defects  [10,  11]  and  interface  incompatibility 
in  multi-phase  composites  [49],  A  Somigliana  dislocation  may  be  idealized  in 
terms  of  a  closed  volume  V  bounded  by  surface  S  located  within  the  interior 
of  a  body  B  of  infinite  size,  as  shown  in  figure  3.  The  surface  S  is  first  cut,  and 
the  two  faces  of  the  cut  undergo  a  small  relative  displacement  b(x),  where  x  is 
the  position  of  the  point  on  the  surface  S.  When  b(x)  produces  separation  of 
matter,  the  empty  space  is  then  filled.  On  the  other  hand,  when  b(x)  produces  an 
overlap  of  material,  we  consider  the  extra  material  to  be  subsequently  removed. 
The  faces  of  S  are  cemented  together,  such  that  any  further  displacements  are 
continuous  across  S.  The  body  containing  the  Somigliana  defect  is  now  in  a 
state  of  internal  stress  (see  also  [50]). 

If  b(x)  =  constant,  the  Somigliana  dislocation  is  equivalent  to  the  usual  transla¬ 
tional  dislocation  of  Volterra  [41]  and  Burgers  [80].  If  b(x)  =  co  x  r(x),  where  co  is 
the  spatially  constant  Frank  vector  [81]  and  r  is  the  disclination  radius,  the 
Somigliana  defect  is  equivalent  to  the  disclination  of  Volterra  [41]  and  Frank 
[81].  To  model  a  point  defect,  we  let  S'  be  a  small  sphere  with  a  distribution  of 
b(x)  over  the  surface.  Letting  the  radius  of  S  shrink  to  zero,  while  at  the  same  time 
prescribing  b(x)  such  that  the  displacement  at  a  fixed  distance  from  the  sphere 
remains  finite,  provides  a  mechanical  representation  of  a  point  singularity  [11].  For 
the  particular  case  of  b(x)  constant  and  directed  radially  outward  from  S ,  we 
obtain  a  model  of  a  spherical  interstitial  or  substitutional  atom,  while  a  vacancy 
is  modeled  for  b(x)  directed  inward.  By  considering  a  uniform  distribution  of 
somewhat  larger  hollow  spheres,  we  can  model  porosity  associated  with  microvoids 
in  crystals. 

We  now  make  the  association  of  the  general  micromorphic  framework  of 
(21)  (28)  with  the  physics  of  crystal  defects.  This  is  accomplished  by  considering 
three  limiting  cases:  traditional  (geometrically-necessary)  dislocation  theory, 
characterized  by  Qcha  =  0;  disclination  theory,  characterized  by  Qcha  —  Qc[ba] ;  and 
a  model  for  isotropically-distributed  point  defects,  characterized  by  Qcha  —  QcCj;a. 
Superposition  of  the  three  cases  then  enables  us  to  consider  Somigliana  dislocations 
of  arbitrary  character. 
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5.1.  Dislocations 

The  description  afforded  by  our  framework  is  restricted  to  translational  dislocations 
if  we  prescribe  Qcha  —  0,  for  which  (22)  reduces  to  the  crystal  connection  [19]: 


r..d  _ 

ch  —  1  cb 


pLa 

x  .a 


t?L—  \a 
*.btc 


_  jfLa  i fL — lc* 

—  .a,  c”  .b 


(29) 


The  torsion  tensor  of  (29)  is  then  written  as 


p..d  _  p..d  _  pLd  pL' 


1  a 

[b,c]  > 


(30) 


which  vanishes  according  to  (35)  when  the  lattice  deformation  is  holonomic,  i.e. 

ffb  =  0  Ffcl?  =  0  **  F'a-]u  =  3axf  (31) 


Figure  4  provides  a  two-dimensional  visual  interpretation.  On  the  left  side  of 
figure  4,  we  inscribe  two  orthogonal  fields  of  constant  contravariant  ‘lattice  slip 
vectors’  on  B  —  denoted  by  u  and  v  —  that  form  a  rectilinear  grid.  Such  a  repre¬ 
sentation  of  a  ‘perfect’  lattice  in  B  is  in  agreement  with  classical  crystal  plasticity 
theory  [82],  for  example,  wherein  Fp  is  assumed  to  leave  the  slip  directions  unper¬ 
turbed.  Notice  that  the  lattice  may  have  slipped  with  respect  to  its  position  relative  to 
the  reference  configuration  B0 ,  consistent  with  figure  2(b)  and  indicated  by  the 
dotted  unit  cells  on  the  left  side  of  figure  4  (for  clarity,  these  slip  steps  are  not 
shown  on  the  right  side  of  figure  4,  though  they  do  exist  in  configuration  B).  Since 
the  lattice  slip  vectors  are  spatially  constant  over  the  relaxed  crystal,  we  have 


u(d)  =  u  (B),  v(A)  =  v(B),  (32) 

where  A  and  B  are  two  neighboring  points  in  the  relaxed  configuration  that  are 
separated  by  a  small  distance  dx  =  x(A)  —  x(  B)  when  mapped  onto  the  deformed 
coordinates  of  the  crystal.  Focusing  for  the  moment  on  the  field  u,  we  multiply  both 
sides  of  the  first  of  (32)  by  F L(A),  the  value  of  the  lattice  deformation  at  the  location 
corresponding  to  point  A ,  to  obtain 

F!f(A)if(A)  =  FLi(A)ud(B).  (33) 


configuration  B  configuration  B 


Figure  4.  Visualization  of  the  crystal  connection. 
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Expanding  F L(A)  with  a  first-order  (i.e.  linear)  approximation  in  x,  we  have 


FLa(A)  »  FlZ(B)  +  FL£h(B)dxh,  (34) 

which  is  then  substituted  into  equation  (33)  to  yield 

ua(A)  =  u“(B)  +  FL£h(B)FL-]a(B)dxbuc(B),  (35) 

where  we  have  defined  the  deformed  slip  vectors  as  push-forwards  via 

ua(A)  =  FL£(A)u*(A),  ua(B)  =  FL£{B)if(B).  (36) 

Using  the  identity  (F^aF^~la)  h  =  0,  we  can  rewrite  equation  (35)  as 

0  -  ua(A)  -  ua(B)  +  FL£(B)FL-f(B)dxbuc(B).  (37) 

Making  the  linear  approximation 

ua(A)  -  u\B)  «  dU^  (xh(A)  -  xh(B))  ,  (38) 

Y 

=dxb 


and  dropping  the  notation  signaling  localization  to  point  B ,  equation  (37)  then 
becomes 


0  =  m"  hdxh  +  FL£’FLc£blddxhuc,  (39) 

which  obtains  the  structure  of  a  null  covariant  derivative,  i.e.  parallel  transport: 

Vhua  =  u%  +  rfcuc  =  0,  (40) 

since  dxh  is  now  arbitrary.  Since  the  Christoffel  symbols  corresponding  to  V  in  (40) 
are  equivalent  to  those  given  in  equation  (29),  parallel  transport  of  the  slip  vectors 
with  respect  to  the  crystal  connection  physically  corresponds  to  construction  of  the 
deformed  lattice  from  a  rectilinear  grid  in  relaxed  configuration  B.  Analogously,  the 
field  v  in  figure  4  satisfies 


V/,v°  =  vah  +  rhacvc  =  0,  (41) 

i.e.  the  vector  field  v  =  FLv  e  TB  constructed  by  pushing  forward  from  the  uniform 
reference  grid  v  e  TB  obeys  the  rule  of  parallel  transport  with  respect  to  the  con¬ 
nection  r.  It  is  clear  why  the  crystal  connection  is  said  to  possess  the  property  of 
‘distant  parallelism’  or  ‘teleparallelism’. 

The  closure  failure  B  of  the  area  a  (with  local  differential  one-form 
ncda  =  ecahdxa  A  dxh,  where  nc  denotes  a  unit  vector  normal  to  a)  shown  in 
figure  4 — i.e.  the  area  of  the  parallelogram  enclosed  by  u(v4),  u(5),  v(A),  and  v(C), 
with  C  a  third  location  on  the  deformed  crystal — is  equivalent  to  the  area  integral 
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of  the  field  of  local  Burgers  vectors  b(x)  [53]: 

B“=  ba  da  —  aahnhda=l-  ffdxh  a  dxc.  (42) 

a  a  ^  a 

where  2a ,ab  =  shdc  fd“  is  the  geometrically  necessary  dislocation  density  tensor 
referred  to  B.  In  terms  of  straight  defect  lines,  the  dislocation  density  tensor  is 
written  [30,  74] 

a  =  Pj V  ®  (43) 

j 

with  pJ,  bJ,  and  the  line  length  per  unit  current  volume.  Burgers  vector,  and 
tangent  line  in  the  spatial  configuration  for  dislocation  population  j.  To  deduce 
(42)  directly  from  the  right  side  of  figure  4,  we  begin  with  the  rather  obvious  vector 


addition  relation 

B“  =  u\A)  -  ua(B)  +  va(C)  -  va(A).  (44) 

Also  from  the  current  configuration  lattice  depicted  in  figure  4  we  have 

(«)  (v) 

dxh  =  x\B)  -  x\A)  =  vb(A),  dxb  =  xb(C)  -  xb(A)  =  ub(A).  (45) 
Upon  making  the  linear  approximations 

(«)  (>0 

u“(B)  -  ua(A)  »  uah(A)dxb ,  va(C)  -  va(A)  »  vah(A)dxb ,  (46) 

and  invoking  equations  (40)  and  (41),  equation  (44)  becomes 

(«)  (v)  _  («)  _  (v) 

B"  =  -u%(A)dxb  +vab(A)dxb  -  rfcdxb  if  {A)  -  rfcdxb  vc(A),  (47) 

which  upon  usage  of  (45)  can  be  written  as 

Ba  =  2f[£]  vb(A)uc(A)  =  f£(v(A)u(A)fc] .  (48) 

Making  the  identification 

(v(A)u(A)fc]=  i  (vb(A)uc(A)  -  vc(A)ub(A)j  =  i  dxb  A  dxc,  (49) 

equation  (48)  is  then  rewritten  as 

B "  =  1  Tbacdxb  A  dxc,  (50) 


which  is  the  local  form  of  (42).  Integrability  conditions  for  the  crystal  connection 
coefficients  of  equation  (29)  have  a  natural  geometric  interpretation,  as  demon¬ 
strated  in  previous  works  by  Schouten  [55]  and  Le  and  Stumpf  [4],  Partial 
coordinate  differentiation  of  (29)  yields 


/  jfL,—  lot  _  t?L—\(x 

\ /  -a  1  be  )  ?  d  r  .c,  bd  ? 


(51) 
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the  left  side  of  which  is  expanded  to  read 


pt-  1  c 
*  .a,  d 


r..a  |  17  b — la  p..fl  17  b — lc 

be  +  f  .a  1  bc,d  ~  ^  .a 


( rfehl  +  hld )• 


(52) 


Since  the  order  of  partial  differentiation  on  the  right  side  of  (51)  is  arbitrary, 
we  then  have 


0=  -2  F 


■L-\d 
c,  [bd\ 


=  F 


L—\a 


(rfcj 


r£d  +  rferdc 


rdaerhec)  =  f 


L—\ci 


R;, 


ebeb 


(53) 


where  R'^d  are  the  components  of  the  Riemann-Christoffel  curvature  tensor  derived 
from  using  (10).  Upon  multiplication  of  (53)  through  by  FL  we  arrive  at  the 
conclusion  that  Rcgd  —  0  are  conditions  ensuring  satisfaction  of  (29).  From  satisfac¬ 
tion  of  (53),  the  crystal  connection  is  said  to  be  ‘integrable’  [55].  Since  the  crystal 
connection  is  metric  with  respect  to  CL  (from  equation  (24)  with  Q£.  =  0),  and  its 
curvature  tensor  vanishes,  the  set  {B,t,  CL}  constitutes  a  metric,  non-Riemannian 
space.  We  also  have 


^  fiJL  t~i.  .ci  ✓ — 'L  1  /^b  s^b 

^bsabl  ( cd)  —  ^ cb,d  1  ^ db,c  ^ dc,b> 


(54) 


easily  verified  by  direct  calculation.  If  we  now  further  impose  that  the  torsion 
T;%  -  0,  such  that  =  F^,  meaning  that  the  connection  is  symmetric  (i.e.  non- 
Cartan  space),  and  physically  that  geometrically  necessary  dislocations  are  absent, 
then  (54)  defines  the  Levi-Civita  connection  coefficients  ( cf.  (8))  on  the  current 
configuration  with  Riemannian  metric  tensor  CL,  and  the  set  {B,  t,  CL}  constitutes 
a  Euclidean  space. 


5.2.  Disclinations 

Volterra  [41]  introduced  six  fundamental  types  of  defects  in  elastic  bodies:  three 
types  of  translational  displacement  discontinuities,  known  as  edge  and  screw 
dislocations,  and  three  types  of  rotational  incompatibilities,  later  termed  disclina¬ 
tions  by  Frank  [81]  and  further  classified  as  either  wedge  or  twist  disclinations 
(see  figure  5).  Disclination  theory  has  been  applied  to  numerous  problems  of  interest. 
These  include  descriptions  of  micropolar  rotations  in  liquid  crystals  [81],  rotational 
defect  substructures  and  commensurate  strain  hardening  in  metal  forming  processes 
[83],  grain  boundary  structures  [84],  deformation  twins  [69],  and  polycrystalline 
triple  junctions  [85].  Disclinations  have  also  been  recognized  as  characteristic 
defects  in  polymers  [86]  and  nanocrystals  [87].  Molecular  dynamics  simulations 
incorporating  disclination  concepts  [88,  89]  have  been  undertaken  to  characterize 
energy  distributions  over  a  range  of  misorientations. 

Continuum  theories  of  distributed  dislocations  and  disclinations  can  be  found  in 
the  geometrically-oriented  papers  of  Anthony  [43],  Eringen  and  Claus  [73],  Lardner 
[44],  Minagawa  [45,  46],  Amari  [47]  and  De  Wit  [48].  Precise  mathematical  descrip¬ 
tions  and/or  elastic  solutions  are  also  included  in  the  texts  of  Nabarro  [42],  Lardner 
[74],  Mura  [90],  Maugin  [1]  and  Zubov  [91].  An  extensive  review  of  disclination 
theory  focusing  upon  defect  kinetics  and  contributions  to  plastic  strain  hardening 
was  provided  by  Seefeldt  [78],  who  suggested  that  partial  disclination  dipoles  be  used 
to  describe  misoriented  subgranular  interfaces  manifesting  upon  grain  refinement 
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Figure  5.  Volterra’s  defects:  reference  cylinder  with  defect  line  §  and  cut  surface  S  (a);  edge 
dislocations  (b,  c)  and  screw  dislocation  (d)  with  Burgers  vector  b;  twist  disclinations  (e,  f )  and 
wedge  disclination  (g)  with  Frank  vector  to. 


in  advanced  stages  of  plastic  deformation.  Li  [84,  92]  demonstrated  for  a  volume 
of  isotropic  linear-elastic  solid  material  the  equivalence  between  the  strain  energy 
of  a  wedge  disclination  dipole  and  a  finite  wall  of  edge  dislocations.  Kroner  and 
Lagoudas  [93]  developed  a  gauge  theory  of  disclinations  intended  to  model  micro¬ 
rotations  in  liquid  crystals.  Lazar  and  Maugin  [94]  solved  for  the  stress  field  of 
a  wedge  disclination  using  higher  gradient  elasticity  theory. 

The  coefficients  rcl  in  the  dislocation-disclination  version  of  our  framework  are 
given  by  equation  (22),  with  the  covariant  components  of  Q  anti-symmetric  [45,  46]: 


Qcba  —  Qcb^cla  Qcab  Qc\ba\' 


(55) 


The  micro-rotation  kinematic  variable  Q  describes  spatial  gradients  of  lattice  rota¬ 
tion  within  the  crystalline  volume  element  arising  from  distributed  disclinations. 
Consider  parallel  transport  of  a  lattice  director  da,  over  the  small  distance  dx  to  a 
new  orientation  da,  conducted  with  respect  to  the  covariant  derivative  (21).  For 
illustrative  purposes,  we  temporarily  restrict  attention  the  special  case  when 
FL  »  1  and  spatial  gradients  of  the  lattice  deformation  are  small.  Applying  the 
summation  convention  over  covariant  Cartesian  coordinates,  we  may  then  write 

d a  —  T  da IfdXfr  —  da  T  Qboc^^b^c  —  ( $ac  T  Q.b[ac\^--^b)^c  —  (56) 


where  Qbac  —  Qb[ac ]  from  equation  (55),  and  the  rotation  matrix  d>ac  satisfies 
t|>7  %  <f>  1  and  det(O)  ~  1  for  small  magnitudes  of  Qbacdxh,  i.e.  small  relative  lattice 
rotations  such  as  those  occurring  across  low-angle  grain  or  subgrain  boundaries. 
We  see  from  (56)  that  Qhac  acts  as  an  effective  ‘gradient’  of  rotation  in  the  spatial 
direction  xb .  However,  since  Q  is  not  required  to  satisfy  the  following  compatibility 
conditions,  it  is  not  a  true  spatial  gradient,  i.e. 


Qb\ac\,d  Qd[ac],b  ^  ^  Qb\ac\  ^\oc].  h' 


(57) 


4002 


J.  D.  Clayton  et  al. 


where  d[ac]  is  a  skew  matrix  that  can  exist  only  when  the  conditions  Qb[ac],  d  —  Qd[aci  b 
are  met. 

We  now  return  to  the  general  situation  wherein  components  of  VL  may  be 
arbitrarily  large.  From  (24),  the  connection  is  metric,  since  the  covariant  derivative 
of  CL  vanishes: 


VcC’ab  =  ~2Qc(ah)  =  0.  (58) 

Components  of  the  torsion  tensor  T  of  the  connection  T  are  given  by 

T£  =  f£  +  2Q$l,  (59) 

where  7):)'  is  the  torsion  of  the  crystal  connection  discussed  in  section  5.1  [16,  19,  21]. 
The  components  of  the  Riemann-Christoffel  curvature  tensor  Rb"d  are,  from  (28) 
and  (55), 


R[ab]cd  —  R(ab)cd  —  2.V[cQd\(ab)  +  TcdQe(ah)  —  0. 


(60) 


Thus  we  see  that  the  curvature  is  limited  to  nine  independent  components,  i.e. 

Rabcd  —  R[ab][cd\  • 

Consider  a  Burgers  circuit  c  in  the  current  configuration,  enclosing  area  a 
comprised  of  oriented  elements  n da.  A  total  Burgers  vector  accounting  for 
incompatibility  induced  by  torsion  and  curvature  may  be  written  as  [40,  44,  45] 


B“  =  - s' 


dbc 


'  be 


'  R-ecb^ 


’^7 tdda  = 


=  Bn 


B‘i 


(61) 


where  2Btj  =  sdhc \  Tbcndda  describes  the  closure  failure  of  c  and 
2 B“r  =  sbdc\  Recabxendda  measures  the  change  in  direction  of  position  vector  xe 
upon  parallel  transport  about  c  with  respect  to  the  connection  f.  Equation  (61) 
reduces  to  (42)  when  disclinations  are  absent  (i.e.  when  T£  =  7};",  Q)'“  —  0,  and 
R'bfc  —  0)  .  We  can  re-write  (61)  in  terms  of  the  second  rank  geometrically  necessary 
dislocation  tensor  a  and  second  rank  geometrically  necessary  disclination  tensor  0, 
each  referred  to  the  spatial  configuration  B: 


B“ 


(aad  +  CL-Ufefgb6sdxb)ndda , 

a  ^  ' 


(62) 


where 


2aad  =  edhcfhac,  4&d  =  f?ba  edce  Rahce.  (63) 

Figure  6  illustrates  the  total  Burgers  vector  B  =  Bx  +  BR  introduced  in  equation  (62) 
in  terms  of  parallel  transport  of  a  lattice  director  vector  if,  about  an  incompatibility 
circuit  c  [40].  When  viewed  from  the  standpoint  of  a  Volterra  process  [41],  one 
can  imagine  the  body  in  figure  6  to  consist  of  the  superposition  of  a  single 
edge  dislocation  (producing  the  incompatibility  BT)  and  a  single  wedge  disclination 
(producing  the  incompatibility  BR). 
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Figure  6.  Burgers  vector  with  contributions  from  dislocations  (Bx)  and  disclinations  (BR). 


Defect  tensors  a  and  0  contain  information  enabling  us  to  fully  reconstruct 
T  and  R: 


T'cb  —  T[cb ]  —  scbdO‘ad ,  R-abcd  —  R[ab][cd\  —  £baeecd^ ■  (64) 

The  tensors  of  (63)  are  also  related  to  the  summed  contributions  of  discrete  crystal 
defects  as 


a  =  A'  ®  \J,  0  =  ®  C* .  (65) 

j  k 

with  p',  b',  and  If  the  net  scalar  dislocation  line  density.  Burgers  vector,  and  unit 
tangent  line,  respectively,  for  dislocation  population  j,  and  with  tjk,  cok,  and  if  the  net 
scalar  disclination  line  density,  Frank  vector,  and  unit  tangent  line,  respectively,  for 
disclination  population  k. 

Notice  that  a  and  6  do  not  account  for  curved  defect  segments  and  combinations 
of  defect  lines  that  do  not  contribute  to  a  net  Burgers  vector.  This  is  evident  when 
(65)  is  written  in  terms  of  positively-  and  negatively-signed  defect  populations  [29]: 

a  =  (p|  -  P-)K  ®  V,  e  =  X!  {p+  ~  rl-)<a+  ®^k’  (66) 

j  k 

where  b+  =  — bi  and  wk+  —  —call.  It  is  emphasized  that  while  the  net  defect  densities 
of  (63),  (65)  and  (66)  represent  incompatibility,  the  toted  mobile  defect  densities 
contribute  to  the  rate  of  plastic  deformation  (to  FP),  the  latter  which  may  be  nonzero 
even  in  the  absence  of  net  defect  densities,  as  would  occur  under  conditions 
of  homogeneous  plastic  flow  (e.g.  when  pJ+  —  pL).  We  use  the  qualifiers  ‘net’, 
‘geometrically  necessary’,  and  ‘non-redundant’  interchangeably  to  describe  defect 
densities  producing  nonzero  incompatibility. 

Linearized  compatibility  equations  for  the  defect  density  tensors  of  (63) 
follow  from  identities  of  Bianchi  and  Schouten  (c/.  [55]),  expressed  in  Cartesian 
coordinates  as 

T[bc,di  =  Am  (cL<*)tb=  e“bcOhf  hw,e]  =  0  0abb  =  0.  (67) 

In  a  small-strain  formulation  that  employs  additive  elastoplastic  strains  and  rotation 
gradients,  De  Wit  [48]  inferred  from  equations  analogous  to  (67)  that  disclinations 
may  act  as  sources/sinks  for  dislocations,  and  that  disclination  lines  cannot  end 
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abruptly  within  a  crystal.  When  disclinations  are  absent,  the  first  of  (67)  depicts  a 
divergence-free  dislocation  density  tensor,  meaning  that  dislocations  cannot  start 
or  end  within  the  crystal. 


5.3.  Point  defects 

We  focus  attention  now  upon  a  geometric  representation  of  uniformly-distributed, 
isotropic  point  defects.  The  same  connection  T  listed  in  equation  (22)  is  used,  but 
the  form  of  Q  is  restricted  as  follows  [45]: 


r-£  =  FfF 


La  t?L —  lot 


Qcl  =  r'd  +  QA  Qcb  =  QX ,  Qcba  -  QcQ 


ba’ 


(68) 


From  (68),  the  contribution  from  point  defects  to  the  derivative  S7cdadxc  is  given 
by  Q'Cadxc  —  ShaQcdxc,  a  pure  stretch.  The  connection  is  non-metric  since,  by  (24) 
and  (68), 


VcC^  = 


■2  Qc(ab)  —  _2  QcCab, 


(69) 


meaning  physically  that  vacancies,  interstitials,  or  micro-voids  alter  the  counting 
of  atomic  steps  in  the  crystal  [102].  The  torsion  becomes,  upon  substitution  of  (68) 
into  (25), 


T'..a  T~>..a  jn ..a  rri..a  ,  /y  pa  c'a 

ch  =  1  cb  ~  1  be  —  1  cb  +  <Jcdb  ~  <Jb°c, 


(70) 


thereby  including  a  contribution  to  the  non-redundant  dislocation  density 
2 aad  =  edhc  Tf  from  the  crystal  connection  and  a  contribution  from  micro- 
morphic  contraction/expansion  in  the  form  Qh8ac  —  Qc8ah.  From  equation  (27), 
the  Riemann-Christoffel  curvature  tensor  becomes  simply 


hed  =  (2%Qd\  +  TfdQe)sab  =  2Qld,c]8ah-(2F[cedlQe  -  TcedQ))8ab  =  2 8ahQ[dc].  (71) 

V 

=0 


Upon  parallel  transport  around  a  current  configuration  circuit  c,  a  lattice  director 
vector  da  changes  by  an  amount  Ad„  as  follows  in  the  presence  of  a  non-vanishing 
curvature  tensor  R'bfd  [44,  55,  56]: 


Adfl  =  d^  -  dfl  =  - 


kicAbdxc  A  dxd  =  -  d 


Q[c,d\dxC  A  dxd 


(72) 


where  we  have  used  (71).  From  the  last  of  (72),  we  see  that  the  vector  dfl  is  unchanged 
in  direction  upon  completion  of  the  circuit  c;  its  components  merely  undergo  stretch 
or  reduction. 

Recall  that  the  micromorphic  variable  Q  represents  additional  degrees-of- 
freedom  at  the  microscale  for  capturing  the  effects  of  defects  on  spatial  gradients 
of  the  lattice  directors.  Accompanying  Q  at  the  macroscale  is  the  contribution  of 
these  defects  to  F  of  (17).  The  contribution  of  isotropic  defects  may  be  written  in 
terms  of  our  F1  term,  neglecting  here  the  contributions  of  anisotropic  residual  lattice 
deformation: 


Fi  =  (u«5)ga®^, 


(73) 
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where  the  scalar  u  is  given  by 

u  =  (1  -  0)-1/3,  (74) 

where  <p  will  be  defined  shortly.  Consider  differential  volume  elements  dV  e  B  and 
dV  e  B.  The  relationship  between  these  volume  elements,  upon  appealing  to  (73)  and 
(74),  is  then 


dV  —  J'dV  =  det  F'<J  det  g/  det  g  dV  =  dV( 1  —  <p)  1 ,  (75) 

upon  assuming  coincident  coordinate  frames  on  B  and  B  (and  thus  equivalent 
metric  tensors  g  and  g).  Equation  (75)  then  provides  the  definition  for  the  scalar  f: 

</>=  (dV-dV)/dV.  (76) 

When  microvoids  or  interstitials  are  considered,  cp  is  the  fraction  of  defects  in  the 
intermediate  configuration,  per  intermediate  configuration  volume  [95]: 

<pdV=  (N+V^  +  dV)  -dV=N+V+,  (77) 

where  A and  V*  are  the  number  of  defects  and  the  volume  of  each  defect,  respec¬ 
tively,  within  the  local  volume  element  dV,  in  configuration  B.  For  an  isotropic 
distribution  of  vacancies,  there  is  a  net  decrease  in  volume  due  to  F1,  and  we  have  [96] 

</>  dV  —  {—AT  V~  +  dV)  —  dV  —  —ATV~,  (78) 

with  N~  the  number  of  vacancies  within  volume  element  dV ,  each  inducing 
reduction  V~ . 


6.  Summary 

Our  theory  is  founded  upon  two  major  assumptions,  the  first  being  a  three-term 
decomposition  of  the  average  deformation  gradient  for  a  crystal  element  (see  also 
equation  (17)): 

F=  r^Fp,  (79) 

f£ 

with  the  F1  term,  non-standard  in  the  usual  theories  from  the  literature  (cf  [97]), 
accounting  for  the  presence  of  defects  that  affect  the  average  lattice  arrangement  and 
internal  residual  stresses  within  the  crystalline  volume  element  (see  figure  2).  Also 
in  (79)  are  the  elastic  deformation  Fe  representing  both  the  recoverable  lattice  stretch 
associated  with  the  average  applied  stress  acting  on  the  element  and  rigid-body 
rotations  of  the  lattice,  as  well  as  the  plastic  deformation  Fp  accounting  for  the 
partition  of  fluxes  of  mobile  defects  that  leaves  the  lattice  unperturbed.  The  total 
lattice  deformation  is  written  as  FL.  The  second  major  assumption  is  an  additive 
decomposition  of  a  linear  connection  describing  spatial  gradients  of  the  slip 
directions  and  lattice  director  vectors  between  neighboring  crystalline  elements 
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(i.e.  the  microscopic  description).  Christoffel  symbols  of  this  connection,  introduced 
in  equation  (22),  may  be  written  as 

f?b=F%I^+CL-ladQ(ibdi+  +CL~lad{Qc(hd)  -  QcCLhd),  (80) 

crystal  micro-rotation  micro- expansion  micro-strain 

connection  ( disclinations )  ( isotropic  point  defects )  ( Somigliana  dislocations ) 

(GNDs) 

where  the  first  term  on  the  right  side  describes  gradients  of  the  director  vectors  due  to 
first-order  gradients  in  the  average  lattice  deformation  tensor  FL,  following  Bilby 
et  al.  [19],  as  discussed  in  section  5.1.  The  micromorphic  variable  Q  (c/.  [45,  48]) 
participates  in  the  remaining  three  terms  on  the  right  of  (80):  a  micro-rotation 
associated  with  disclinations,  as  discussed  in  section  5.2,  an  isotropic  micromorphic 
expansion  associated  with  point  defects,  discussed  in  section  5.3,  and  a  general 
micromorphic  strain  that  may  be  used  to  represent  arbitrary  lattice  director 
deformations  when  superposed  with  the  other  terms  in  (80).  Dislocation  and 
disclination  density  tensors  then  follow  from  the  torsion  and  curvature,  respectively, 
of  the  connection  (80),  the  latter  vanishing  when  the  variable  Q  vanishes.  Our  theory 
is  unique  from  others  in  the  literature  via  its  combination  of  the  three-term  deforma¬ 
tion  decomposition  (79)  with  the  additive  connection  for  geometrically-necessary 
defects  (80).  Note  that  it  is  possible  to  describe  defect  kinematics  (e.g.  construct 
the  dislocation  and  disclination  density  tensors)  without  decomposing  the  lattice 
deformation  VL  into  distinct  recoverable  and  residual  terms.  For  example, 
Minagawa  [45,  46]  invoked  a  connection  similar  to  (80),  but  used  only  a  single  elastic 
lattice  deformation  map,  amounting  in  essence  to  a  two-term  deformation  gradient 
decomposition.  However,  we  find  our  three-tenn  decomposition  more  realistic  than 
the  usual  two-term  decomposition  from  the  viewpoint  of  developing  constitutive 
relations,  with  only  the  recoverable  portion  Fe  of  the  lattice  deformation  VL  acting 
as  a  thermodynamic  force  conjugate  to  the  externally  applied  stress  [67]. 

In  summary,  we  present  table  1,  which  depicts  the  classifications  of  the  affine 
geometries  corresponding  to  crystal  defects  in  finite  elastoplasticity,  following 
terminology  of  our  section  2  and  motivated  by  Steinmann  [6].  Figure  7  is  a  two- 
dimensional  idealization  of  a  crystalline  element  in  the  spatial  configuration 
containing  the  lattice  defects  listed  in  table  1.  We  remark  that  in  many  figures  in 
this  paper,  simple  cubic  lattices  have  been  used  to  illustrate  key  concepts,  as  this 
crystal  structure  is  easiest  to  visualize.  Our  framework  has  been  developed  with  the 
primary  intention  of  describing  engineering  metals  that  deform  by  dislocation  glide, 
most  typically  exhibiting  a  Bravais  lattice  structure  and  most  often  realized  in  a 
face-centered-cubic,  body-centered-cubic,  or  hexagonal  close-packed  arrangement. 
The  concepts  forwarded  herein  are  most  naturally  applied  to  cubic  lattices,  wherein 
the  director  vectors  of  (18)  may  be  assigned  parallel  to  the  primitive  translation 
vectors  of  the  lattice.  The  framework  is  valid  for  crystal  structures  of  lower  symmetry 
(e.g.  hexagonal)  so  long  as  the  Cauchy-Born  rule  ( cf.  [103])  applies  for  the  elastic 
deformation  of  the  interatomic  bond  vectors.  In  such  cases,  while  the  directors  of 
relations  (18)  assigned  to  each  crystal  volume  element  are  not  parallel  to  the  physical 
edges  of  the  unit  cell  of  the  lattice,  they  are  assumed  to  fully  characterize  the  stretch 
and  rotation  of  the  interatomic  vectors  comprising  the  primitive  cell.  On  the  other 
hand,  for  more  exotic  structures  wherein  the  recoverable  deformation  is  non-uniform 
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Table  1.  Geometric  classifications,  kinematic  quantities,  and  corresponding  crystal  defects. 


Geometry 

Torsion 

Curvature 

Metric 

Defects 

Euclidean 

T  =  0 

R  =  0 

VCL  =  0  (yes) 

Statistically  stored 

(B  admits  holonomic 

dislocations 

coordinates) 

T^O 

R  =  0 

(SSDs)*  figure  7a 

Cartan  and  non- 

VCL  =  0  (yes) 

Geometrically  necessary 

Riemannian 

dislocations  (GNDs) 

(e.g.  f  =  f) 

T  =  0 

R^O 

figure  7b 

Riemannian  and 

VCL  =  0  (yes) 

Disclinations 

symmetric 

T^O 

R^O 

Metric,  Cartan, 

VCL  =  0  (yes) 

GNDs  and  disclinations 

and  Riemannian 

T  fO 

R^O 

figure  7c 

Non-metric,  Cartan, 

VCL  f  0  (no) 

All  types  of  defects 

and  Riemannian 

(e.g.  point  and  line 
defects)  figure  7d 

*  Admissible  defects  for  all  geometries  in  table  1.  However,  all  dislocations  are  geometrically 
necessary  when  considered  in  isolation;  collections  of  SSDs  produce  no  net  Burgers  vector 
and  no  net  torsion  tensor. 


5,  r,  T=  0 


B,  r,  0 


Figure  7.  Current  configuration  lattices  accompanying  table  1:  (a)  Euclidean  space, 
(b)  Cartan  space,  (c)  Cartan  and  Riemannian  space,  (d)  general  non-metric  space. 


within  the  unit  cell,  additional  degrees  of  freedom  not  included  in  our  theory 
are  required — see  for  example  the  couple  stress-based  treatment  of  diamond 
cubic-structured  silicon  by  Garikipati  [98]. 


6.  Conclusions 

A  comprehensive  and  rigorous  theoretical  framework  for  the  finite  deforma¬ 
tion  kinematics  of  lattice  defects  in  crystalline  materials  has  been  presented. 


4008 


J.  D.  Clayton  et  al. 


The  framework  includes  a  three-term  multiplicative  decomposition  of  the  deforma¬ 
tion  gradient,  with  terms  accounting  for  recoverable  elasticity,  residual  lattice 
deformation  due  to  defect  populations,  and  plastic  deformation  resulting  from 
defect  fluxes.  Also  invoked  is  an  additional  micromorphic  variable  representing 
additional  degrees-of-freedom  associated  with  lattice  defects  other  than  conventional 
translational  dislocations,  such  as  disclinations,  point  defects,  and  most 
generally,  Somigliana  dislocations.  Our  treatment  has  been  purely  kinematic  in 
nature;  accompanying  thermodynamics  and  kinetics  remain  to  be  addressed  in 
future  work. 
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